Two-Fluid Atmosphere for Relativistic Stars 

E.N. Glass* and J.P. Krisch 

Department of Physics, University of Michigan, Ann Arbor, Michgan 48109 

(December 21, 1998) 



"q^" ■ We have extended the Vaidya radiating metric to include both a radiation 

Q\ ■ fluid and a string fluid. This paper expands our brief introduction to exten- 

t— i ; 

£H . sions of the Schwarzschild vacuum which appeared in 1998 Phys. Rev. D 57 

R5945. Assuming diffusive transport for the string fluid, we find new analytic 
solutions of Einstein's field equations. 

> 

m . PACS numbers: 04.20.Jb, 04.40.Dg, 97.60.-s 

o 
o 

On 
On 

"o 
cr 

> 

X 



*Permanent address: Physics Department, University of Windsor, Ontario N9B 3P4, Canada 



1 



I. INTRODUCTION 



Strings have become a very important ingredient in many physical theories. They may 
have been present in the early universe and played a role in the seeding of density inhomo- 
geneities [0. The idea of strings is fundamental to superstring theories f2j. String fluids 
have been suggested as dark matter candidates. The lensing properties of cosmic strings 
and string systems have been discussed || and the apparent relationship between counting 
string states and the entropy of the Schwarzschild horizon @], |5[], || suggests an association 
of strings with black holes. Recently Glass and Krisch (?]] have pointed out that allowing the 
Schwarzschild mass parameter to be a function of radial position creates an atmosphere with 
a string fluid stress-energy around a static, spherically symmetric, object. If the mass is also 
a function of retarded time a Vaidya radiation fluid is present in addition to the string fluid. 
Since the metric has one arbitrary function, m(u,r), invariantly defined by the sectional 
curvature of the u, r two-surfaces, a given mass distribution determines the stress-energy. 

The string fluid stress-energy is a macroscopic, statistically averaged, description of a mi- 
croscopic distribution of Planck length string bits. The string fluid lies on a two-dimensional 
timelike world-sheet (u, r) and its stress-energy is parametrized by radial derivatives of the 
mass function. The radial stress, p r < 0, is the expected string tension. Since the string 
fluid is constrained to stay in the (it, r) plane, the Planck scale string bits have only radial 
motions. Interactions of the string bits along the radial directions are modeled by the macro- 
scopic string tension p r and create, on average, a long string stretching radially away from 
the core or horizon of the mass distribution. The Glass-Krisch atmosphere allows transverse 
stresses. These stresses, if present, can be modeled by the presence of a pressureless dust 
fluid. The transverse pressures arise from the different velocities of the dust and string fluids. 
The most general atmosphere has three components: a Vaidya radiation fluid, a string fluid 
with radial tension, and a dust fluid. This can be interpreted as a two-fluid atmosphere with 
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a transverse pressure. We also discuss other, less general, possibilities. 

This atmosphere can model a variety of physical situations at different distance scales. 
It can describe the atmosphere around a black hole with a distance scale of multiples of 
Schwarzschild radii. It could also describe a globular cluster with a dark matter component 
and a distance scale of the order of parsecs. 

The contracted Bianchi identities are satisfied for arbitrary m(u, r) and so the mode 
of mass propagation is a modeling choice from the many possible propagation equations of 
classical and quantum physics. Each choice allows the generation of new analytic solutions to 
the field equations. In this work we develop new solutions by considering diffusive transport 
of the string fluid elements. With this choice we find that the string fluid diffuses inward as 
the Vaidya photons carry energy outward. The overall effect is to slow down the time scale 
of the Vaidya energy loss. 

In the next section we present our extension of Vaidya' s metric [§] as a further extension 
of the Schwarzschild vacuum solution. The arbitrary mass function, which appears in the 
metric, is identified as the sectional curvature of nested two-spheres. The string fluid and 
an associated static two-fluid model are described in the third and fourth sections. We are 
able to provide a two-fluid kinematic interpretation of the transverse stresses in this model. 
In the fifth part of the paper diffusive transport is developed and analytic solutions for the 
energy density are found. In the sixth section analytic mass solutions are presented and 
interpreted. The last two sections contain the analysis of the horizon structure of some of 
our solutions and a discussion of the general results. 

In this work Greek indices range over (0,1,2,3) = (u, r, </?). Our sign conventions are 
2-4i/;[ Q /3] = A^R^ ua/3 , and R^ v = R a liua . Overdots abbreviate d/du, and primes abbreviate 
8/ dr. Overhead carets denote unit vectors. We use units where G — c — 1. Einstein's field 
equations are G^ v = StcT^, and the metric signature is (+,-,-,-). 
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II. EXTENDING THE SCHWARZSCHILD VACUUM 



The spacetime metric covering the region exterior to a spherical star is given by 

ds 2 = Adu 2 + Idudr - r 2 (dtf 2 + sin 2 ^ 2 ) (1) 

where A = 1 — 2m(u,r)/r. Initially m(u,r) = mo provides the vacuum Schwarzschild 
solution in the region r > 2rriQ. At later times that region admits a two-fluid description 
of radial strings and outward flowing short-wavelength photons (sometimes called a "null 
fluid"). Metric ([!]) is spherically symmetric and given in retarded time coordinate u. With 
the use of a Newman-Penrose null tetrad the Einstein tensor is computed from ([!]) and given 
by 

Gfxv = -2®u(lfinv + n^l u + m^fh^ + fh^m^) (2) 
-2$ 22 ^^ - 6A^. 

Here the null tetrad components of the Ricci tensor are 

$ n = (2m' - rm")/(4r 2 ), (3a) 
$ 22 = -m/r 2 , (3b) 
A = R/24 = (rm" + 2m') /(12r 2 ). (3c) 

The only non-zero component of the Weyl tensor is 

M/ 2 = -m/r 3 + (4m' - rm")/(6r 2 ). (4) 

The metric is Petrov type D with and principal null geodesic vectors 

l^dx^ = du, (5a) 

n^dx fl = (A/2)du + dr, (5b) 

m^/ = -(r/V2)(cM + i sinfidtp), (5c) 

where 



l K u = (A' /2)l fl l v - {l/r)(m^m v + m^m v ), (6a) 
n^-u = -(A' /2)n lx l l/ + (A/2r)(m lx m l/ + m^m u ), (6b) 
m Kl/ = (A/2r)l tl m u - {\jr)n jX m v + (coti?/ \[2r) {m^m u - m^rh u ). (6c) 

In order to clearly see the two-fluid description we introduce a timelike unit velocity vector 
■0 M and three unit spacelike vectors f M , y3 M such that 

The unit vectors are defined by 

v^dx? = A l ' 2 du + A~ 1/2 dr, = A~ 1/2 d u , (7a) 

TpdaP = A-^dr, = A~ 1/2 d u - A l/2 d r , (7b) 

d^dx 11 = rd$, = -r -1 ^, (7c) 

ft dx ,i = rsm&dip, (p^d^ = -(rsini?) -1 ^. (7d) 

is hyper surface-orthogonal, i.e. v\p, ;v v a ] — 0, with h^ v the first fundamental form of the 
hypersurface. Since v^dx^ = f(u,r)dt, the components of show explicitly that -0 M lies 
along t = const time lines: 

h tlv dx p 'dx v = (g^ u — v il v v )dx^dx v (8) 
= -A- X dr 2 - r 2 {d$ 2 + sin 2 i?V)- 

The kinematics of the flow are described by 

v\ = a»v v + o* v - (6/3) (f"f „ + W„ + <p*<p v ), (9) 



where 



= [m/r + A9 r (m/r)]A- 3/2 f M , (10a) 
< = (6/3) (-2^ + W„ + $><p v ), (10b) 
= (m/r)A~ 3/2 . (10c) 



The Einstein tensor can now be written as a two-fluid system: 

G M „ = (2m/r 2 )y F - (2m'/r 2 )(v^v u - f^ v ) (11) 
+K/r)(^i + ^). 

Spherical symmetry allows the function m(u, r) to be identified as the mass within two- 
surfaces of constant u and r, and invariantly defined from the sectional curvature [28] of 
those surfaces: 

- 2m/r 3 = Rvp^tffrp. (12) 



III. STRING FLUID 

The string bivector is defined by 

E " = £ WTSfi> (^ = (0,1) or (2,3). 

Spherical symmetry demands that the averaged string bivector will have a world-sheet in 
either the (u, r) or (p) plane. The condition that the world-sheets are timelike, i.e. 
7 := lE^E^j, < 0, implies that only the S nr component is non-zero. It is useful to write 
Tj^ u in terms of unit vectors 

TT = r»v v - v^r v . (13) 



It is now clear that S^E^ = v^v v — f^r v '. We follow Letelier ||, [10] and write a string 
energy-momentum tensor by analogy with one for a perfect fluid 

T^ id = p Ufi u u -ph^, 

where h^ u = — u^u vi h^ u u u = 0. The string energy-momentum is given by 

T string = p( _ T )l/^ af^ _ p±H (M) 
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where H^ v = — E^E^, H^JiyP = 0. Although here 7 = —1, we have kept 7 explicit in 
(14]) and written E Ml/ := (— / y)~ l / 2 Y l flu to show how E^ is made invariant to reparameteri- 
zations of the world-sheets || . 

Einstein's field equations allow the matter portion of G^ u in Eq.([TT|) to be identified as 
a string fluid: 

T, u = T™ d + T s ;j m °, (15) 
= i>lplu + pv^Vu + p r r^r u + p±(-& ll '& v + 0p0 u )- 

Thus 

Atci/j = —rh/r 2 , (16a) 
4np = —4iip r = m'/r 2 , (16b) 
87rp_i_ = — m" 'jr. (16c) 

Since the contracted Bianchi identities are satisfied for arbitrary m(u,r), it follows that the 
equations of motion = are identically satisfied for the components of T Mi , given in 
Eq.©. 

The components of the contracted Bianchi identities which do not vanish because of 
spherical symmetry but rather because of the explicit form A = 1 — 2m(u,r)/r are 

= -V„[(2$ n + i?/4)r] - 

and 



n^!i = -V,[(2$ n + R/A)n v + $ 22 Z"] - G^n 



IV. STATIC FLUID MODELS 

Static models are useful in developing insights about time dependent fluids. We consider 
two static models that are equivalent to the general stress-energy tensor described in Eq.(p~5|) 
when it is time-independent. The first is a static isotropic string fluid, and the second 
provides a two-fluid interpretation of (|15D and an interpretation of the transverse stress. 
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A. Isotropic string fluid 



Consider the stress-energy tensor in Eq.(|TjJ) with static mass function m(r) and with 
Pr = p±: 

T ™v = ~Pr(VnVv - r p r v - - (17) 

This is clearly an isotropic cloud of strings with equation of state p + p r = 0. The mass is 
determined by the constraint p r = pj_ or 



m" nt 



(18) 



2r r 

The solution of Eq. (|l8D , in the more general time-dependent case with null radiation fluid, 



is 



m(u, r) = r 3 ci(u) + c 2 (tt). (19) 

For the static case we have c\ and c 2 constant. As can be seen from energy-momentum 
Eq.fll7|) the isotropic string cloud is in an Einstein spacetime. 

B. Static two-fluid model 

One can use two different four- velocities, and w^, to write a two- fluid model: 

T 2flmd = p2l b^w v + P\UyU v + p r f„f v . (20) 

The fluid with p 2 is dust and the other is a fluid with a radial stress, (m = and there 



is no Vaidya radiation fluid.) Letelier [|Tj| has described a procedure for casting two-fluid 
stress-energies into the form of an anisotropic fluid. His method is not adapted for the string 
fluid equation of state. As a variation of his method, we transform the fluid velocities to 
create two un-normalized vectors: timelike and spacelike. 

Vfh = yfpl cosa Uf, + ^fp~ 2 sina (21) 
v /p 2 " = -yfpi sina ii^ + Jp~ 2 cosa w^. (22) 
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This transformation obeys 

Pxu^Uv + P2W p w v = piV/jVu + p 2 Y fl Y u . (23) 

Since is spacelike, Eq.(^) is valid only for w M ^ and a 7^ n{ji/2). The parameter a 
is defined by V^Y 1 * = 0: 

cot2a = — ^ 1 ^ 2 — . 
2y/Pip2 « M W 

The stress-energy T^} utd can be written as 

T tran S / = p ^ + p ^ + (24) 

with 

P := ^ Q V Q Pl , p ± := —Y a Y a p 2 . (25) 
The transformation in Eq.(|23|) constrains the new density and transverse pressure to obey 

pi + p2 = P~P±- (26) 

T transf hag ^ game form ag rpstring ^ Eq.(^). 

If the two fluids have the same velocity, = w^, then the two-fluid stress-energy becomes 
a single fluid with a radial stress. 

T 2fluid ^ rplfluid = {pi + + pr ^ K (27) 

The density of this string fluid is just the sum of the two densities. In the case w M = the 
transverse pressure is zero. The transverse pressure reflects the different velocities of the 
fluids in the two- fluid model. When the transverse pressure is zero then m = c\t + ci- 

Using a multi-fluid model, one can also describe the radial stress as reflecting a velocity 
difference between two dusts so that the complete stress-energy tensor ( [3~5|) could be de- 
scribed by a three-fluid dust model with all of the dusts moving with different velocities. 
Which of the interpretations discussed in this section is most likely depends on the actual 
physical situation being modeled. 
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V. DIFFUSIVE TRANSPORT 



Diffusion has been a seminal process for the development of our understanding of many 
modern systems. Since the work of Einstein \ Wjj , Smoluchowski [17], and Langevin [18], the 
ideas implicit in diffusion have found new application in areas as diverse as the behavior 
of stock option values to cosmic strings. Vilenkin jOJ has introduced diffusion into the 
description of cosmic strings by characterizing string evolution as the formation of Brownian 
trajectories. This has also been discussed by Bennett P0fl. Diffusion is playing an important 



part in the growth of our understanding of fluctuations in quantum gravity |21| , f22| and in 



the very early life of our universe. For example, Watabiki |23| has used the classical diffusion 
equation to characterize diffusion times in fractal quantum gravity [^4] , [ 25 1 . Diffusion may 
also play a role in eternal inflation models where inflation field fluctuations can be modeled 



as random walks J26J. The path integral technique, developed by Norbert Weiner |]27J to 
describe diffusive processes, has become an essential part of the modern view of quantum 
mechanics. 

String collisions, unlike point particle collisions, do not occur at a single spacetime point 
(interaction vertex). The two-surface picture of strings requires the collision (interaction) to 
be a curve on a world-sheet. Observers in different Lorentz frames will see the interaction 
occurring at different points along the curve. Statistically, the coarse-grained picture of 
phase space for a set of collisions hides the lack of Lorentz invariance of a single collision. 
We assume string diffusion is like point particle diffusion in that the number density diffuses 
from higher numbers to lower according to 



d u n = V V 2 n 



(28) 



where V 2 = r~ 2 (d / dr)r 2 (d / dr) and T> is the positive coefficient of self-diffusion (which we 
henceforth take to be constant). Classical transport theory derives the diffusion equation 
by starting with Fick's law 



J(n) = -VVn 



(29) 
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where V is a purely spatial gradient. Then 4-current conservation JLy^ = 0, where 

J (n) d V = (n, J(n)) (30) 

= nd u — V(dn/dr)d r , 

yields the diffusion equation (|28|) . We label the 4-current J( n ) to indicate n diffusion but we 
could have also written Ji p \ since the string number density and string fluid density must 
be related by p = M s n where M s is the constant mass of the string species. M s must be 
a multiple of the Planck mass since it is only over Planck length scales that point particles 
resolve into strings. 

By rewriting the components (|16a|) and (|16b| ) as m = — 4nr 2 ifj and m' = 4irr 2 p, we 



can write the integrability condition for m as 

p + r- 2 d r (r 2 ilj) = 0. (31) 
If we compare the diffusion equation Q2E| ) (n replaced by p) 

p = V r- 2 d r (r 2 dp/dr) (32) 

with p in Eq . (|3lD we obtain 

m = AtcV r 2 dp/dr. (33) 

Thus solving the diffusion equation for p and then integrating those solutions to obtain m 
provides exact Einstein solutions which can be interpreted as either anisotropic fluids or 
diffusing string fluids. 

There are some analytic solutions of Eq.(^2|): 



P = Po + h/r, (34) 

P = Po + (k 2 /6)r 2 + k 2 Vu, (35) 

p = p + k 3 (Vuy 3/2 exp[-r 2 /{AVu)}, (36) 

P = Po + (&6/ r )exp(— klVu)[sm(k4r) + k 5 cos(k i r)}. (37) 
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Solutions ([34] ) and (|36D appear in J7|. In addition to the solutions above, we have obtained 
a number of analytic solutions using well known similarity techniques. Those solutions will 
be presented elsewhere. The physical behavior of the density solutions provides a variety of 
atmospheric models. 

Solution fl5%D describes an atmosphere with a simple drop off in radius and no time 
dependence. k\ must be positive to avoid negative densities, po is the density at spatial 
infinity. 

The density described by equation (|35| ) has two very different behaviors. For &2 > 0, the 
density is not physically realistic. It increases with radius and grows with time. However, 
for k 2 < 0, the density decreases with radius and has a zero indicating a bounded string 
atmosphere. Since we are working with a string fluid whose equation of state is p = —p r , 
the boundary also has zero radial pressure. The position of the boundary moves inward as 
T>u increases; the extent of the atmospheric shell decreases with time. 

The third density solution, Eq.(^), requires k 3 > for positive densities. The solution 
begins with a high central density which falls off with radius. As time progresses the density 
decreases to the constant value p , the density at spatial infinity. 

The fourth density solution, Eq. Q57D , models a complete array of atmospheric behaviors. 
For example, the parameter choices k^/po — 10, k^ = k$ = 1, describe an atmosphere that 
starts with a string boundary at T>u = 1 and as time progresses becomes an unbounded string 
cloud. Other parameter choices model atmospheric shells which are always unbounded. 
Parameters k 5 and k e must be positive for positive densities. 

VI. DIFFUSIVE MASS SOLUTIONS 

A. Analytic solutions 

Upon integrating m' = 4irr 2 p and fn = 4nV r 2 dp/dr we obtain , in the same order as 
the densities above, the following 



12 



m(u, r) = m + (4vr/3)r 3 po + lixk^r 2 - 2Vu), (38) 

m(u, r) = m + {An/3)r 3 p + (4tt /3)r 3 k 2 {Vu + r 2 /10), (39) 

m(u, r) = m + (47r/3)r 3 po + I6irk 3 [-r]exp(-r] 2 ) + ( v / 7r/2)erf(?7)], (40) 

m(u, r) = m + (47r/3)r 3 p + (4:rckQ/kl)exp(-klVu)B , (41) 

where i] := r (42>u) _1//2 , and -B = sin^r) — A;4rcos(/c4r)+/c5[cos(/c4r) + /c4rsin(fc 4 r)]. Solutions 
(|38l) and (fiTf ) appear in 0. 

If parameters fci, £;2, and k§ are zero then all the mass solutions above become the 
static solution m + (47r/3)r 3 p - This is the mass for the isotropic string fluid described by 
Eq.©. 



B. Interpreting the mass 

Diffusive transport has a conserved four-current for the density 

J^d^ = pd u - V{dp/dr)d r 

which can be written in terms of the null tetrad 

r = pn» + [{A/2)p - V(dp/dr)]F (42) 

When J! 1 = is integrated over a four-volume, Stokes theorem casts the integral onto the 
bounding three-surfaces: 

/ J^^9d A x = j> J^V^gdS, (43) 

Ri 8R 4 

where we integrate over u = const null three-surfaces A/2 and Af± with dS^ = l^drdfidip, and 
r = const timelike three-surfaces E 2 and Si with dS^ = (n^ — ^l^dudddip. The four-volume 
i?4 can be pictured as a tube surrounding a timelike cylinder containing the central source. 
An orthogonal cross-section of i? 4 would have sides A/i, A/2, Ei, E 2 , forming a rhomboid 
with Ei bounding i? 4 away from the central source. Thus 
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/ J^l^-^gdrdMip + J J^(n^ - ^l^y/^gdud&dtp = 
A/b-A/i S2-S1 
with specific form 

J psf^gdrdtidip - J V(dp/dr) y f=gdud'dd<p = 0. (44) 
A/i— M S2-S1 
To first understand the string mass, we examine a static string fluid. Metric (|l|) includes 

static string fluids when m(u,r) is restricted to m(r). We can use the density solution 

Eq.(|34"D with D = so that there is no diffusion. For p = p + fci/r, we have 



4tt / (p r 2 + k x r)dr = 0. (45) 
A/2 -A/i 

Thus the string mass is static with m str i n g(r) = J = J 



A/i M 



2 



m st ring(r) = m Q + 4vr(^p ^ 3 + ^i^ 2 ) (46) 

on null surface A/i and at a later time on J\f 2 . This is the mass in Eq. (|38D with T> = 0. 

Now consider the case with time dependence where the mass is m(u, r) and there is a 
Vaidya fluid of short photons. Again using p = po + ki/r, we find 

4vr J (p r 2 + k x r)dr + 4vr J k{Ddu = (47) 
A/b-A/i S2-S1 
This yields the total mass of strings and null fluid in Eq.(B3) where 



m(u,r) = J - J = J - J (48) 

A/2 Ei A/i £2 

We use our knowledge of the static case to identify the string mass as that part of the total 
mass integrated over the null three-surface A/", and identify the flux through the timelike 
surface £ as resulting from both the energy carried by the Vaidya photons and the diffusing 
strings. Photons enter R4 from the central source through Si while string bits diffuse through 
Ex toward the source, with the opposite happening at E 2 . There the photons leave i? 4 while 
string bits enter. Thus m(u, r) = m string + rrifi ux where 

m s tring = An(^p r 3 + ^k x r 2 ) 
rrifiux = -AnkiVu 

Of course m str ing need not be static. The density examples with p(u, r) will have m string (u, r). 
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VII. HORIZONS 



The topological two-spheres ip) nested in an r = const surface at time u have outgoing 
null geodesic normal and incoming null geodesic normal n M . The spheres become trapped 
surfaces when both l^r,^ and n M r )A1 are positive (negative for +2 signature). The marginally 
trapped surface is the outer boundary of all trapped surfaces at time u, and the apparent 
horizon is the time history of the marginally trapped surface. 

Here Pr, M = 1. The fluid atmosphere has ?i/V, M = —A/2, negative until m(u,r) > r/2. 
At that time the mass has become compact enough to trap light. 

We will analyze mass expression ([38]) for trapped surfaces. With Lq := 2mo, L^ 3 : = 
(47rpo)/(3mo), L^ 2 '■— 27iki/mo, and 1/Tq := AnkiD / m®, we search for zeros of the polyno- 
mial 

at time u = const. The four parameters (To, Lq, L\, L?) are positive. Suppose there is a root 
at r = aL , a > 0. Then, at u — 0, Eq.(^9|) becomes 



a 3 (L /L 1 ) 3 + a 2 (L /L 2 ) 2 + 1 - a = 



or 



a 2 (2m ) 2 [a(87rp /3) + 27rA;i/m ] + 1 - a = 0. (50) 

If fci and po are zero then a = 1 and we find the Schwarzschild horizon at r = 2niQ. For 
non-zero parameters Eq. ( |50|) has no real roots, which is consistent with the outgoing short- 
wavelength photons. 

At time u = and for some short time after 1 — u/Tq is positive. After more time Am, 
1 — u/T becomes negative and then there is only one sign change in Eq. (|49|) . Descartes' 
rule of signs tells us that, for u > Au, there is at most one real root of Eq. (|49|) and the fluid 
will have a trapped surface at r = 2m(u, r). 
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VIII. CONCLUSION 



That t> M is hypersurface orthogonal implies v^dx^ 1 = f(u,r)dt. If metric (p]) were the 
Schwarzschild metric then A = 1 — 2rrio/r, u = t — r — 2moln(r — 2mo), and v^dx^ = A~ x l 2 dt. 
Next in simplicity is the Vaidya metric, where A of metric ([I]) would be A = 1 — 2m(u)/r. 
The Vaidya metric is not static and one cannot coordinate transform back to Schwarzschild. 
With the Vaidya metric r — 2m(u) is a spacelike hypersurface lying outside the local null 
p9| , unlike the Schwarzschild hypersurface r — 2mo which is null. 
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The traditional view of the Vaidya metric places it outside a spherical star which is 
losing mass via Vaidya's short-wavelength photons. Vacuum Schwarzschild geometry is 
joined smoothly to Vaidya at its radiative boundary. 

The system described here with m(u,r) continues to have f) M hypersurface orthogonal, 
so the t in v^dx^ = f(u, r)dt labels spacelike hypersurfaces with energy-momentum given 
by Eq . (|T5|) above. Because of the implicit equation of state, p + p r = 0, we interpreted the 
fluid as an "atmosphere" of open strings which form a string fluid. This view has also been 
discussed by 't Hooft || and Maldacena |J. This fluid was taken as a source of energy- 
momentum. As a modeling choice, we took the strings to interact diffusively with energy 
carried away by short-wavelength photons. It was necessary to use the approximation of 
short- wavelength photons since (2m/r 2 )I /1 l 1/ is not an exact solution of Maxwell's equations. 

Quantum field theory requires the quantum vacuum to be Lorentz invariant [ 30(1 , which 



constrains the energy-momentum tensor on a Minkowski background to have the form = 
pr]^ u , and the energy- density to transform as 

, = p + p(v 2 /c 2 ) 

p i _ ( v 2 /c 2) 

under change of inertial frame. The string equation of state p r = —p satisfies the required 
transformation property and so the vacuum outside a relativistic star could indeed include 
a string atmosphere. 

This work necessarily forms an incomplete picture of the evolution of the astrophysical 
system we model here. Strings may exist at the Planck length scale, with a large number 
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of them possibly providing a macroscopic classical string, and forming a visible atmosphere 
until they lie within a trapped surface. 
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